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A theoretical analysis was performed focusing on the heat transfer associated with 
the laminar oscillating flow in a tube. This situation finds applications in the Stirling 
engines or regenerative-type refrigerators where the working fluid in the heat 
exchanger undergoes an oscillatory motion. Such an oscillating flow conceivably 
entails the thermally developing region, not only because the swept length of working 
fluid is roughly equal to or longer than the characteristic length of the heat exchanger, 
but also because the wall temperature changes abruptly along the longitudinal direc¬ 
tion. For simulation of the practical heat exchanger composed of cooler and heater, 
two types of thermal boundary conditions are taken into account; either wall tempera¬ 
ture or wall heat flux has a square-wave distribution. It is found that the thermally 
developing length increases in proportion to the oscillation frequency at slow oscil¬ 
lation but eventually approaches an asymptotic value at high frequency. The local 
average Nusselt number in the developing region is observed to be inversely pro¬ 
portional to the square root of the distance measured from the thermal discontinuity. 
Out of the thermally developing region, the local Nusselt number is determined only 
by the oscillation frequency regardless of axial position. © 1998 Elsevier Science Ltd. 
All rights reserved 
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Nomenclature 

Nu, 

local Nusselt number 



Nu a 

Nusselt number in case of a sinusoidal wall 

c 

n,k 

function of A defined in Equation (7) 


temperature distribution 

C /m d m 

complex coefficient 

Pr 

Prandtl number 

D 

diameter of pipe [m] 

q w 

dimensionless wall heat flux defined in 

^m,n 

function of 77 and A defined in Equation (11) 


Equation (14) 

8k 

temperature profile function defined in 

R 

radius of pipe [m] 


Equation (6) 

r 

radial coordinate [m] 

In 

first kind of the modified Bessel function of 

Re 

Reynolds number 


order n 

T 

temperature [°C] 

i 

imaginary unit number, V - 1 

T 

x a 

amplitude of a sinusoidal wall temperature 

7,i 

first kind of the Bessel function of order n 


distribution [°C] 

L d 

thermally developing length [m] 

t 

time [s] 

4 

swept length [m] 

u 

velocity [m/s] 

L w 

period of the wall temperature distribution 

u m 

mean velocity [m/s] 


[m] 

X 

dimensionless axial coordinate, 2irz/L w 

Nu 0 

Nusselt number in case of a sinusoidal wall 
temperature distribution as A—>0 

z 

axial coordinate [m]. 
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CO 

angular velocity [rad/s] 






absolute value 






Real[ ] 

real component of a complex number 

Greek letters 

0 



J 

L 

* 

P 

modified Womersley number, R^IcoIk 

section averaged quantity, (j) 

- = 2 

fydy 

V 

dimensionless radial coordinate, r/R 



0 

0 _ 

K 

thermal diffusivity [m 2 /s] 



1 

11 

f 


A 

ratio of the swept length to the period of an 


time averaged quantity, / = 

/dr /277 


arbitrary wall temperature distribution 



J 

< 



A 

ratio of the swept length to the period of a 


locally averaged quantity over the axial 


sinusoidal wall temperature distribution 



* 



A c 

e 

critical swept length ratio 
dimensionless temperature, T/T a 


direction, f z - 

fdz/z 



0 w 

dimensionless wall temperature 


0 



Ad 

dimensionless temperature difference 

1 

radial gradient at 17 = 1,/ = 

df/drj( 1) 

T 

dimensionless time, cot 

sign( ) 

sign of a real number; sign(v) = + 1 for v > 



0; sign(v) = - 

1 for v < 0. 




Introduction 

The flow pattern in the heat exchanger of a regenerative- 
type cycle machine such as a Stirling engine/refrigerator 
and a Vuilleumier refrigerator displays the oscillating fea¬ 
tures in which the working fluid reverses its flow direction 
periodically in time. Furthermore, the relevant fluid dynam¬ 
ics and heat transfer phenomena in these system mostly 
occur in the hydrodynamic ally and thermally developing 
regime because the length of the heat exchanger is usually 
shorter than the swept length of the working fluid and the 
wall temperature changes abruptly along the longitudinal 
direction. Besides, the pressure inside the system oscillates 
as well, due to the variation of the internal volume that the 
fluid occupies in the machine and/or due to the heat transfer 
between the fluid and the heat exchanger. Owing to the 
complexity involved in the phenomena, the heat transfer in 
this situation has not been analyzed comprehensively con¬ 
sidering all these effects and still remained unanswered. 

In the early days of research, most of the investigations 
treated the heat transfer phenomena rather simple way in 
their analysis; an incompressible oscillating flow with lin¬ 
eally approximated wall temperature distributions 1-3 . They 
were successful to reveal the basic features of the heat 
transfer involved in the oscillating flow with the above 
assumption. More recently, some researchers have included 
the compressible effect of the flow field in the analysis 4 - 5 . 
In doing so, they approached the problem with linearization 
of the governing equations which adopt only the first har¬ 
monics and discard all other higher harmonic terms. While 
the linearization of the governing equations made the prob¬ 
lem relatively simple for analytical investigation, this 
approach could still be readily applicable to the realistic 
situations where the wall temperature distribution can be 
linearly approximated or the swept length of the oscillating 
flow is much shorter than the characteristic length of the 
heat exchanger. Thus, for example, these investigations 
have contributed enormously to the better understanding of 
the heat and work flow especially in the thermoacoustic 
devices. 

However in the heat exchanger of a Stirling cycle 
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machine, the linearization approach can not be applied 
since the swept length of the working fluid is in the same 
order of magnitude as the length of the heat exchanger and 
furthermore, the wall temperature is not linear anymore but 
rather changes abruptly along the longitudinal direction 6 . 
Therefore, the major portion of the heat exchanger pertains 
to the hydrodynamically and thermally developing region. 
Consequently, it is necessary to examine the heat transfer 
within the thermally developing region to properly explain 
the heat transfer phenomena in the Stirling cycle machine. 

Several articles have been published about the 
developing region of the oscillating flow but are confined 
to the fluid flow itself 7 - 8 . The studies over the attendant heat 
transfer in the developing region have been few. 

Though the flow condition is not exactly the same as the 
oscillating flow, it is worthy to note the theoretical work of 
Siegel and Perlmutter 9 on the heat transfer in the thermally 
developing region of a pulsating flow between parallel 
plates. Assuming uniform velocity over the channel cross 
section, they investigated the explicit dependence of the 
overall heat transfer on the pulsating frequency. Recently, 
Kim et a/. 10 earned out numerical simulations based on a 
more realistic model and explored the principal mechanism 
of the heat transfer by the pulsating flow. However, the 
above-stated works cannot deal with the situations where 
there exist flow reversals due to the assumptions applied in 
the analysis. Hence, it is not certain that their results could 
be extended to the extreme case in which a steady compo¬ 
nent of the pulsating velocity is zero, i.e. the oscillating 
flow. 

Therefore, our purpose in this paper is to examine theor¬ 
etically a heretofore unrecognized aspect of the heat trans¬ 
fer process in the thermally developing region by an oscil¬ 
lating flow in a circular pipe. In this interest, the general 
solution has been found firstly for an arbitrary wall bound¬ 
ary condition by superposing the temperature distributions 
developed for the sinusoidal wall temperature distri¬ 
butions 11 . Then, applying the general solution to the situ¬ 
ations with the thermal discontinuity along the pipe, the 
temperature field is obtained in the thermally developing 
region and the dependency of the thermally developing 
length on the oscillation frequency is examined. Also, the 
attendant heat transfer characteristics in the thermally 
developing region are investigated. 

Although this work deals with an incompressible oscil- 
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lating flow, we do believe that the current research would 
contribute to the better understanding of the heat transfer 
in a Stirling cycle machine since the careful examination 
has been exercised to accommodate the characteristics of 
the thermally developing region which has never been 
investigated yet. 


I 0 (P^ikr)) 

IoipJik) 


Q*( A) = (-!)* +1 4(2A)/„_*(2A). 


( 6 ) 


(7) 


Analysis 

Sinusoidal wall temperature distribution 

In order to set the scene for the rest of this article, the 
fundamental case with a sinusoidal wall temperature distri¬ 
bution is adopted from Lee et a/. 11 and summarized. By 
applying the assumptions of constant thermal properties, 
negligible axial conduction, and uniform velocity over the 
cross section with pure sinusoidal oscillation (i.e. u = 
L s co cos ((ot)/2), the governing equation for this case is 
derived as 


_ 7 dO dO 13/ d0\ 

P 2 2 + 2 ^ A cos t 3- = - 5- 77 — , 

OT OX 7) 07) \ 07)I 


do 

0t) 


t, = 0 = 0, 0 T, = 1 = 0w = sinx 


( 1 ) 


where the dimensionless valuables are 


Arbitrary wall boundary condition 

In this section, we construct a solution corresponding to a 
more complicated wall temperature distribution from the 
fundamental temperature distribution, Equation (5), with 
the aid of superposition principle. In this solution pro¬ 
cedure, the assumptions introduced to derive the fundamen¬ 
tal temperature distribution, such as a uniform velocity 
assumption, are consistently effective to the general sol¬ 
ution. Though the uniform velocity profile corresponds to 
a limiting case of Pr<5C 1, the heat transfer characteristics in 
this case do not deviate much from those for a moderate 
value of Pr, thereby the solution with uniform velocity 
assumption can successfully predict the heat transfer in real 
situations with 2D velocity profile 11 . 

In what follows, the general solutions are developed for 
two kinds of wall boundary condition, i.e. arbitrary wall 
temperature distribution and arbitrary wall heat flux distri¬ 
bution. 


0 = 


T 

T 

a 


T = COt, X = 277 




and the nondimensional groups are 


Arbitrary wall temperature distribution. The arbi¬ 
trary wall temperature distribution, whenever it is periodic, 
can be expanded in a Fourier series 


|6 =R 



ttL, 
21 ' 


(3) 


0 W = 2 c,„exp(/mx) 

m = — 00 


In the above equations, T a and L w are the amplitude and 
the period of a sinusoidal wall temperature distribution 
respectively and L s is the swept length of the oscillating 
flow defined as 


in which 


1 

277 


77 

( O w exp(imx) dx. 


j 


( 8 ) 


(9) 


4 


1 

2 


277 

* 

\u\ dr. 

V 

0 


(4) 


The two nondimensional numbers, i.e. /3 and A represent 
respectively a ratio of the radius to the thickness of the 
Stokes’ thermal boundary layer and a ratio of the swept 
length to the characteristic length of the wall temperature 
distribution. 

The solution to Equation (1) is 11 


Since the fluid temperature corresponding to each term 
in Equation (8) can be obtained directly from Equation (5), 
the temperature field subject to the thermal boundary con¬ 
dition of Equation (8) can be expressed as 


0 = 0 W + Real 


2 c m exp(imx)2 G m „(7)A)exp(inT) 


L ra = — 00 


n = 0 


( 10 ) 


0 = sin x + Real ^ C 0k (\)g k sin x 

-k= 1 

+ 2 2 L' 2 n,i<k A)<?/fcexp(2n?T) sinx 

n = 1 Ifc = — oo ; ^#0 

OO 

+ i 2 c 2 n - u(A)g*exp((2« - 1 )/ t ) cos 


k = - o°X#0 


where G,„ „ is 


(5) G„( r|,A) - 


y O I J: ( m.\ n= 0 

k = 1 

OO 

( - 1)" 2 C lhk (mA)g k (7)), n> 0 


k = — 0 


( 11 ) 


in which 


in which A is the ratio of the swept length to the period 
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of wall temperature distribution of Equation (8) and is writ¬ 
ten by 


A = 


ttL, 
IT ' 


( 12 ) 


From Equation (10), the time-averaged fluid temperature 
distribution is reduced to 


d = 0 W + Real 


2 tj,A) exp(wu) 


L m = — oo 


(13) 


where the bar denotes a time-averaged value. 

The wall heat flux is obtained by differentiating Equation 
(10) and evaluating it at r) = 1. The instantaneous wall heat 
flux is then written as 

oo oo 

2 C,„exp( imx)2^G' ,„,„exp( inr) 

m = — °° n = 0 -* 

(14) 


de . 

(^w I ^ ~ i Real 


and the time-averaged value is 


q w = Real 


2 c m G',„ j0 exp(/mx) 


(15) 


where G' is defined as 


s-y r _ U^J mn | 

^ m,n l 1 7 = 1* 


(16) 


Arbitrary wall heat flux distribution. When the wall 
heat flux distribution is given as a boundary condition, the 
given heat flux is usually constant in time. However, if the 
heat capacity of the wall is not zero (i.e. the wall thickness 
is finite), the heat flux at the inner wall varies with time in 
a period due to the effect of the wall absorbing heat from 
and rejecting it to the fluid instantaneously. In other words, 
the given heat rate as the boundary condition is the value 
at the outer wall and the heat flux at the inner wall oscillates 
with time, though its time-averaged value is equal to the 
heat flux at the outer wall. 

Therefore, the boundary condition could be given as the 
time-averaged heat flux at the inner wall and may be 
expressed as 


oo 

q w = 2 d m exp(imx), d {) = 0. 

m — — o° 


(17) 


In the above, d 0 , which is the average heat flux over the 
axial direction, must be zero for the temperature field to 
achieve a periodic steady state. 

When the fluid temperature in this case is expressed in 
the form of Equation (10), the time-averaged wall heat flux 
must satisfy Equation (15). Therefore, a relation between 
the Fourier coefficients for the wall heat flux and those for 
the wall temperature can be obtained by comparing Equa¬ 
tions (15) and (17). Substitution of this relation into Equ¬ 
ation (10) yields the temperature field with given wall heat 
flux distribution as follows: 


1 

0 = 0 W + 2 “ —— - 
oo Real[G mJI ] 

m = oo 

oo 

Real d m exp( imx) ^ G,„ „exp( im) 

n = 0 

where the wall temperature d w is 

^ d m 

= 2j p ,, r , — n exp(imx) 
„ Real[G m , 0 ] 


( 18 ) 


(19) 


Thermally developing region 

Using the foregoing results, we can examine how the heat 
transfer in an oscillating pipe flow is influenced by the sud¬ 
den change of the wall boundary condition along the axis. 

As mentioned earlier, the thermally developing region 
occurs due to the abrupt change of the boundary condition. 
However, a stepwise variation along the axis is not spatially 
periodic and therefore we cannot directly utilize the forego¬ 
ing solutions. Under the situation where the swept length 
of the fluid is much shorter than the characteristic length 
of the wall boundary condition, i.e. 

A«1 (20) 

the square-wave distribution is expected to yield almost the 
same features of heat transfer as in the stepwise distri¬ 
bution. This is due to the fact that the neighboring periods 
of distributions have insignificant influence on the heat 
transfer of the region to be investigated. 

In what follows, we consider two square-wave thermal 
forcings (either temperature or heat flux; as shown in Fig¬ 
ure 1 ) and derive the relevant heat transfer characteristics 
by means of the foregoing generalized solutions. The heat 
is supplied or removed at each plateau in Figure 1 and thus 
each plateau can be physically interpreted as a heater or a 
cooler. The heater and cooler can be represented with two 
different types of boundary conditions; first, the heater and 
cooler maintain the constant wall temperatures in the case 
of temperature forcing boundary condition. Second, they 
supply constant wall heat fluxes in the case of heat flux 
forcing boundary condition. Also, it can be assumed that 
there exists an imaginary piston on the center of each part 
to generate the oscillating flow. 



°r?w 

1 










0 

-1 


z 


t 4 






imaginary piston 


Figure 1 Square-wave-form distribution of the wall boundary 
condition 
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Temperature forcing 

In this case, the wall temperature is expressed as a Four¬ 
ier series 

4 °° 1 

0 W = - 2 -Z - - sin (2m - l)x (21) 

77 2 m — 1 

m = 1 

and then the fluid temperature can be obtained from Equ¬ 
ation (10) as 

4 x, 1 

e = 6 W + - 2 j o -, R eal[G 2 ,„- 1,0 sin (2m - l)x 

77 2m — 1 

m = 1 

00 

+ 2 i G 2,n - i,2«exp(2n/T) sin (2m - l)x (22) 

n = 1 

- iG 2m - i,in- iexp((2« - 1)/t) cos (2m - l)x}]. 

Evaluating from Equation (22), the time evolution of the 
fluid temperature within a period of a single oscillation is 
shown as series of 3D plots in Figure 2. The interval of 


time sequence is selected as 0.577. The first frame of a series 
corresponds to the moment when the fluid arrives at the left 
most side within its cyclic motion and is about to reverse its 
flow direction. The longitudinal coordinate is normalized 
by the swept distance L s instead of L w . With this represen¬ 
tation, the temperature distribution near the thermal dis¬ 
continuity is assured to be independent of the period of wall 
boundary condition as long as A«l. In other words, the 
square-wave boundary condition would yield the same heat 
transfer characteristics with the stepwise boundary con¬ 
dition. 

Considering the case of j8 = 1 in Figure 2, it is found 
that the transverse fluid temperature distribution is para¬ 
bolic and the longitudinal distribution is similar to that of 
wall temperature though it varies a little with time. As (3 
increases, however, the longitudinal temperature distri¬ 
bution of the fluid near the center of the pipe starts to devi¬ 
ate from that of the wall temperature and eventually 
becomes a constant profile which is invariant with time and 
reciprocates as if it is a solid body. This is because at large 
/3, the heat transfer from the wall does not penetrate into 
the core part but is limited only within a thin annular space 
near the wall. Meanwhile, the results obtained in this study 



T = -0.571 X = 0.071 T = 0.571 T = 1.071 

Figure 2 Timewise variations of fluid temperature near the wall temperature discontinuity 
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Figure 2 Timewise variations of fluid temperature near the wall temperature discontinuity 


such that the longitudinal temperature distribution of the 
fluid is completely different from that of the wall at large 
/3, cannot be obtained from the analysis with linealized 
equations and show the distinguished results of the present 
work from others. 

The section-time-averaged temperature of the fluid, (6) 
is shown in Figure 3 to clarify the region affected by the 
thermal discontinuity. 

The thermally developing length may be defined as the 
distance at which the temperature difference between the 
wall and the fluid becomes zero. Then, it can be found 
from Figures 2 and 3 that the thermally developing length 
increases with the oscillation frequency but converges to 
the swept length of the fluid if the frequency parameter /3 
becomes approximately greater than 2. The reason why the 
thermally developing length becomes equal to the swept 
length is that the fluid particles which oscillate in the region 
z > L s cannot ever experience the thermal discontinuity 
(note that the axial diffusion is neglected in the present 
study). Consequently, the effect of the change of the wall 
temperature can expand no farther than the swept length. 
The above discussion is relevant to the result of Iguchi et 
al . 8 that, when the oscillation frequency is high enough, the 




Figure 3 Section-time-averaged fluid temperature in the ther¬ 
mally developing region in case of temperature forcing 
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length needed for the oscillating velocity to achieve a fully 
developed profile is also equal to the swept length of the 
fluid. 

Heat flux forcing 

When the wall heat flux is distributed in a square-wave 
form, the fluid temperature can be obtained by Equation 
(18) similarly to the case of temperature forcing. Unlike 
the case of the temperature forcing, however, neither 0 W 
nor (0) is independent on A over the entire domain. This 
precludes the identification of the thermally developing 
region directly from the temperature distribution of the wall 
or the fluid. Fortunately, it is recognized that the time-aver¬ 
aged temperature difference 

<A 0)=0 w -(0) (23) 


Nu mA implies the Nusselt number in case of a sinusoidal 
wall temperature distribution with A = mA and the charac¬ 
teristics of this quantity are well described by Lee et al. 11 . 

Since we are dealing with the case of A«1 [see Equ¬ 
ation (20)], Equation (25) can be transformed into an inte¬ 
gral by letting (2m - 1)A = A and 2A = dA, such that 


m 


4 

77 


sin 



z 

L s 


(ANu a ) dA. 


(26) 


In order to perform the integration and thus to find the 
asymptotic value of (AO) as z/L s — *►<*> , the functional 
relationship of Nu A upon A is needed. With the reference 
to Lee et. al n , it can be found that Nu A is dependent only 
on (3 when A is small and becomes to increase as A 
increases over a certain value, i.e. 


4 y Real[(G 2 „, _ 

77 ^ (2m - l)Real[G' 2 ,„_ 10 ] 

m = 1 


sin (2m - 1 )x 


is independent of the swept length ratio and thus can be 
utilized to define the thermally developing region. The axial 
distribution of (AO) for several values of /3 is drawn in 
Figure_4. It is evident that there exists an asymptotic value 
of (AO) for each /3, thereby implying the existence of a 
fully developed region at large zJL s . Meanwhile, the vari¬ 
ation of (AO) at small z/L s indicates the presence of the 
thermally developing region. 

The asymptotic value, (A0) asymp , can be found by trans¬ 
forming the summation in Equation (23) into an integral 
and examining the asymptotic behavior of the integral as 
z/L s goes infinite. 

Firstly, Equation (23) can be rewritten with the definition 
of Nu mA 


Real[G' m0 ] 

Real[<G w , 0 >] 


such as 


(24) 


Nu a ~ 


Nu 0 , Nu 0 


lim 

Nu a forA«A c 

A—>0 


1.52/3VA forA»A c 


(27) 


in which Nu 0 represents the Nusselt number for the infini¬ 
tesimally small A and is a function of /3 only. And A c 
implies the critical value of A above which the dependency 
of Nu a upon A changes. The parameter A c can be quantified 
as the value of A at the intersecting point of the two extreme 
curves of Nu A , i.e. 


( Nu 0 \ 

\ 1.52/3/ • 


Since Nu 0 is a function of /3 only, A c also depends only 
on j3. Nu 0 and A c are drawn from Lee et al. n and shown 
in Figure 5. It can be found from the figure that there also 
exist two asymptotes for both of Nu 0 and A c as below: 


Nu 0 


\ 6 as 
[V2jS as 


0-K) 

j8— 


(29) 


(A0) 


8 ^ sin(2m-l)x 

77 ,^1 (2m ~ 1 ) Nu (2m - 1)A ’ 


(25) 



V L , 

Figure 4 Section-time-averaged difference between wall and 
fluid temperatures in the thermally developing region in case 
of heat flux forcing 



Figure 5 Variation of Nu 0 and A c with respect to the modified 
Womersley number 
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With the aid of Nu 0 and A c , the Nu A curves are nor¬ 
malized and depicted in Figure 6. This figure reveals that 
all the Nu a /Nu 0 curves for different values of /3 have almost 
same dependency on A/A c regardless of /3. It is also found 
that the normalized curves have two asymptotes: 

Nu a f 1 as A/A c —>0 

—- « \ _ (31) 

Nu 0 [VA/A C as A/A c ^oo 

Using this relation, Equation (26) can be approximated 
as 


m = 


4 1 

77 Nu 0 




A Nu*\ d /A\ 

A c NuJ \A C / 


(32) 



In order to find the asymptotic value of (AO), we examine 
the behavior of the integrals in Equation (32) when z/L s 
goes infinite. When z/L s » 1, Equation (32) can be simpli¬ 
fied as 12 



0 


=^{ si 8 n<z)+o (( A A) 2 )}' 

Remembering that we are dealing with the case of a 
square-wave distribution of the wall heat flux with the con¬ 
dition of A<SC 1, i.e. the stepwise distribution, Equation (33) 
can be rewritten as 


Here, worthy of remark is the appearance of Nu 0 in Equ¬ 
ation (34). This dictates that the heat transfer coefficient in 
the thermally developed region for the heat flux forcing is 
identical to that corresponding to a sinusoidal temperature 
distribution with small A. Such a similarity for the asymp¬ 
totic behavior of (AO) may arise from the dominance of a 
long-wave-length sinusoidal component over the thermally 
developed region. _ 

Figure 4 reveals that the axial position at which (AO) 
achieves its asymptotic value (i.e. the characteristic length 
of the thermally developing region) increases with (3 for 
small /3 and tends to converge to a constant but apparently 
order of L s . Although this situation is very similar to the 
case of temperature forcing, a quantitative determination 
of the thermally developing length requires an additional 
manipulation due to the tricky behavior of (AO) near the 
transition point. 

The thermally developing length 

Our attempt to quantify the thermally developing length 
begins with introducing a new characteristic length in such 
a manner that all the curves in Figure 4 fall onto a single 
curve, i.e. a normalized curve. A clue to the new length 
scale can be obtained from Equation (32). By choosing 
A C z/L s as a new independent variable instead of z/L s , the 
integrals in Equation (32) become independent of /3 and 
thus the temperature difference curves can be normalized. 

Figure 7 shows the reconstructed plot of Figure 4 adopt¬ 
ing A C z/L s and Nu 0 (AO) as new coordinates. It can be found 
from this figure that all the Nu o (A0) curves for each /3 
coincide nearly with each other, though the profiles are a 
little different. It is also found that the curves converge to 
an aymptotic value all together at the vicinity of A C z/L s = 
1 with a maximum deviation of 2%. This reveals that the 
developing region is confined in the region of A C z/L s < 1. 
Therefore, the length of the developing region can be pre¬ 
scribed as 

L A = ]-L S . (35) 

A c 


(AO) 


asymp 




(34) 


Since A c depends on /3, L d is also a function of /3. The 
variation of the thermally developing length with respect 



0.01 0.1 1 10 100 


UK 

Figure 6 Variation of Nu/Nu 0 with respect to A/A c 



Figure 7 Normalized distributions of the temperature differ¬ 
ence in the thermally developing region in case of heat flux 
forcing 
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to /3 can be obtained by Equation (30) or Figure 5 which 
describes the dependency of A c upon /3. 

In conclusion, the asymptotic behavior of the thermally 
developing length L d can be described as 


f-'d'f-'s 


36 P 

1.52 2 

2 


as /3—>0 


as j3—>°o 


(36) 


The variation of L d /L s together with the asymptotic 
curves from Equation (36) is depicted in Figure 8. It is 
interesting to note that the asymptotic curve for /3—>0 is 
practically valid also for moderate values of /3 and the same 
is true for the upper limiting case of j8—►oo. In practical 
applications, the normalized frequency /3 is an order of 10 
for both heaters and coolers of Stirling cycle machines. The 
results in Figure 8 indicate that the thermally developing 
length of the heater and cooler is nearly equal to the swept 
length of the working fluid. By contrast, the normalized 
frequency /3 is typically below 1 for the regenerator. 
Although our analysis does not account for the presence of 
regenerator, it can be inferred that the thermally developing 
length of a regenerator is proportional to (3r. In view of 
practical applications, our analysis implies that the effect 
of the thermally developing region must be considered for 
the heater and cooler but it is negligible in regenerator due 
to a large disparity in the relevant /3 value. 

It is also worthy to compare the thermally developing 
length between the oscillating flow and the unidirectional 
flow. We consider a low frequency limit, i.e. /3—>■(), since 
the oscillating flow can be thought to have features similar 
to the unidirectional steady flow when the oscillation fre¬ 
quency is extremely low. 

It is well known that the thermally developing length for 
the unidirectional laminar flow in a tube (supposing slug 
flow and constant heat flux with 2% deviation) is approxi¬ 
mately L d /D = 0.058Re Pr 13 . If the mean velocity of the 
oscillating flow 
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is used as a velocity scale in LJD = 0.058Re Pr, we have 
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0.232 , 

L d /L s =- (3 2 . (38) 
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The above equation is identical with that for the case of 
/3—^0 in Equation (36) except the proportional constant. 
This fact states that the basic heat transfer characteristics 
of slow oscillating flow are almost the same as that of the 
unidirectional steady flow. Meanwhile, the coefficient in 
Equation (38) is larger by 15% than that in Equation (36), 
which implies that the true length of the thermally 
developing region by slow oscillating flow is smaller by 
15% than that obtained by quasi-steady assumption. This 
discrepancy comes mainly from the difference in the wall 
boundary conditions. In case of oscillating flow, as dis¬ 
cussed earlier, the heat flux at the inner wall oscillates due 
to the effect of thermal capacity of the wall, while the per¬ 
iod-averaged heat flux is kept constant. Therefore, the heat 
transfer characteristics of the oscillating flow cannot be 
simulated exactly by the quasi-steady assumption presum¬ 
ing the wall heat flux kept constant every instant. Parti¬ 
cularly, in case of heat flux forcing, the instantaneous heat 
flux at the inner wall calculated from Equation (14) is not 
constant in time and not even distributed in a square-wave 
form over the axial direction. 


Heat transfer in the thermally developing 
region 

As discussed earlier, the Nusselt number in the fully 
developed region is a constant equal to Nu 0 . But in the 
thermally developing region, the Nusselt number is a func¬ 
tion of the axial position because the temperature difference 
varies along the axis while the wall heat flux is constant. 

To investigate the effect of the thermally developing 
region on the heat transfer, the local average Nusselt num¬ 
ber, Nu z , is defined as 


Nu. = 2 


q w dz / (A0)dz. 


(39) 


The introduction of the integrated average inherits from 
the need of estimating the performance of a heat exchanger 
when the available length is specified. Given the length of 
a heat exchanger by z, the overall average of the Nusselt 
number for the heat exchanger is estimated to be Nu r 
Figure 9 shows the local average Nusselt number defined 
in Equation (39) for several values of /3. The horizontal 
axis is normalized with the thermally developing length, 
while the vertical axis is normalized with the local Nusselt 
number in the thermally fully developed region, i.e. Nu 0 . 
It is noticeable that the curves make a narrow band in the 
figure, which means that the effect of /3 on the overall heat 
transfer does not vary significantly if the thermally 
developing length is selected as a characteristic length 
scale. It can be also seen that the local average Nusselt 
number is approximately proportional to (z/L d ) _1/2 in the 
developing region and approaches Nu 0 as z/L d becomes 
greater than 1. From the design standpoint, for the heat 
exchanger longer than the thermally developing region, the 
heat transfer is enhanced only in the developing region and 
remains saturated out of the developing region. Therefore, 
the augmentation of the heat transfer rate can be achieved 
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Figure 9 Local average Nusselt number in the thermally 
developing region 

by reducing the length of the heat exchanger, at least 
shorter than the thermally developing length. Nevertheless, 
there must be a trade-off between the heat transfer enhance¬ 
ment and the decrease in the heat transfer area. 


Conclusion 

The heat transfer characteristics of an oscillating flow in a 
tube are investigated based on a theoretical model of a heat 
exchanger comprising a heater and a cooler. The present 
model admits the presence of a thermal discontinuity along 
the axis and therefore is more realistic than the previous 
studies dealing with linear temperature distributions. Two 
types of thermal forcings (i.e. either the temperature or the 
heat flux has a square-wave distribution) are considered as 
the tube wall condition. The solutions are obtained by 
superposing the particular solution pertaining to sinusoidal 
temperature distributions. The results are analyzed and 
presented focusing on the thermally developing length and 
the Nusselt number variation. 

For the temperature forcing, the thermal developing 
length approaches an asymptotic value which is equal to 
the swept distance of the working fluid. 


The heat flux forcing better simulates the thermal con¬ 
dition of a heat exchanger and therefore is analyzed in more 
detail. For the heat flux forcing, the thermal developing 
length increases in proportion to the oscillation frequency 
(i.e. /3 2 ) at low frequency. The participating proportional 
constant is deduced from the similarity of the oscillating 
flow to the unidirectional flow at a low frequency limit. 
Flowever, at a high frequency limit, the thermally 
developing length assumes an asymptotic value of an order 
of the swept distance. Within the thermally developing 
region, the local Nu is inversely proportional to the square 
root of the distance from the thermal discontinuity. On the 
other hand, over the fully developed region, the heat trans¬ 
fer rate is almost saturated and is identical to that corre¬ 
sponding to a sinusoidal temperature distribution with an 
infinitesimally small swept length ratio. In light of the heat 
exchanger design, this confirms that the length of the heat 
exchanger needs to be order of the swept distance. 
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